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A  NESTED- DECOMPOSITION  APPROACH  FOR 


SOLVING  STAIRCASE-STRUCTURED  LINEAR  PROGRAMS 

by 

Philip  G.  Abrahamson 

Abstract 

The  algorithm  solves  a  T-period  staircase-structured  linear  program 
by  applying  a  compac.:  basis-inverse  scheme  for  the  Simplex  Method  in  con¬ 
junction  with  a  choice  mechanism  which  uses  the  dual  of  the  Nested  Decom¬ 
position  Principle  of  Manne  and  Ho  to  determine  the  incoming  basic  column. 

A  sequence  of  one-period  problems  is  solved  in  which,  typically,  information 
is  provided  to  period  t  from  previous  and  subsequent  periods  in  the  form 
of  surrogate  columns  and  modified  right-hand  side,  and  surrogate  rows  and 
modified  cost  coefficients,  respectively. 
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I.  Introduction 


This  paper  describes  preliminary  work  on  an  algorithm  for  solving 
staircase-structured  linear  programs.  Such  problems  often  arise  in  the 
modeling  of  phenomena  which  are  naturally  described  as  evolving  over  a 
sequence  of  temporal  or  spatial  "periods".  A  pure  nested  decomposition 
algorithm  transforms  the  problem  into  an  ordered  set  of  smaller  problems, 
one  for  each  period,  which  are  coordinated  only  through  price  and  acti¬ 
vity  communication  between  adjacent  periods.  The  process  of  achieving  an 
optimal  coordination  involves  repeated  solution  of  the  individual  prob¬ 
lems.  Preliminary  experience  has  indicated  that  the  convergence  of  the 
pure  algorithm  may  be  slow.  To  accelerate  this  convergence,  the  algorithm 
is  modified  to  enable  the  individual  problems  to  communicate  in  an  impli¬ 
cit  fashion  whenever  possible.  This  modification  involves  the  genera¬ 
tion  of  surrogate  columns  which  are  passed  to  subsequent  periods,  allow¬ 
ing  these  period  to  parametrically  adjust  solutions  to  earlier  periods. 

A  compact  basis-inverse  scheme  is  used  to  represent  these  parametric 
variations . 

Section  2  states  the  problem  of  interest  and  describes  the  pure  nested 
decomposition  algorithm.  Section  3  outlines  the  modified  approach  and 
discusses  some  details  of  implementation. 


2.  Nested  Decomposition  of  the  Staircase  Structure 
2 . 1  Mie  Staircase  Structure 
The  problem  of  Interest  is 


minimize 


subject  to 


I  ctxt 
t=l 


Alxl 


~Bt-lxt-l  +  Atxt 


xt  2  0 


t  *  2, 
t  =  1, 


T 

T 


where  xf  is  nt  x  1  ,  is  mt  *  n  ,  and  all  other  vectors  and 
matrices  are  of  conformable  dimension.  This  linear  program  is  said  to 
be  staircase-structured  because  the  constraint  matrix  has  the  form: 


Figure  1:  Staircase-Structured  Constraint  Matrix 
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The  inven- 


Activities  in  period  t  are  represented  by  the  matrix  At  . 
tory  provided  by  period  t  to  period  t  +  1  is  described  by  the  matrix 
-B  .  Otherwise,  the  activities  in  period  t  have  no  direct  effect  on 
either  previous  or  subsequent  periods. 

The  dual  of  (P)  is 

T 

(D)  maximize  £  b 

t=l  C  t 

subject  to  irtAt  -  <  c  ,  t  *  1,  T  -  1 

*1*1  1  CT  ' 


2.2  Nested  Decomposition 

Manne  and  Ho  [7]  and  Glassey  [5]  repeatedly  applied  the  decomposition 
principle  of  linear  programming  developed  by  Dantzig  and  Wolfe  [3]  to 
achieve  a  nested  decomposition  of  (P)  .  A  sequence  of  applications, 
modifications,  and  improvements  has  led  to  advanced  implementations  by 
Ho  and  Loute  [6],  who  have  solved  some  large-scale  problems  more  rapidly 
in  this  fashion  than  by  directly  applying  commercial  linear  programming 
to  (P)  . 

Van  Slyke  and  Wets  [8]  describe  an  algorithm  for  solving  (P)  in  the 
case  T  =  2  which  is  equivalent  to  applying  the  decomposition  algorithm 
to  (D)  .  Dantzig  [2]  outlined  an  algorithm  consisting  of  a  nested  de¬ 
composition  of  (D)  .  This  paper  represents  preliminary  work  on  the  de¬ 
velopment  of  a  technique  which  is  based  on  his  approach. 
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The  algorithm  of  Dantzig  [2]  consists  of  the  following  nested  decom¬ 


position: 


maximize  +  U2b2  +  ^3^  +  ^4^  +‘ 


subject  to  tt^A^  “  ^2Bi 


w2^2  _  1T3B2 


1t3A3  *4B3 


+ 

T  T 


ViVi 


►master. 


<  c. 


<  c. 


master. 


V  . 

I  subl 


^T^T-l  —  CT-lj 
"l^T  -  CT  ! 


1  sub. 


■  master^ 
5T-1 


; 1  t 

]  sub_ 


Figure  2:  Nested  Decomposition  of  (D) 


Let  be  multipliers  for  master^  ^  ,  and  let 

bt  =  bt  +  V1V1  • 


(1) 


Then  the  t-th  period  master  problem  has  the  form 


(Dt) 

maximize 

\bt  +  °tft  + 

subject  to 

\At  +  °tFt  +  PtLt  -  Ct 

:  X 

pte  <  1 

:  0 

0t  ,  f*. 


>  0 
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where  (Ft>  f  )  and  (Lt>  £  )  are,  respectively,  the  extreme  ray  and 
extreme  point  proposals  generated  by  period  t  +  1  .  For  t  =  1  ,  b^  =  b^  , 
and  for  t  =  T  ,  only  the  terms  involving  lr^  are  present.  Note  that 
each  proposal  is  a  row  vector. 

The  dual  of  (Dt)  is 


(pt) 


minimize 


z  =  c  x  +  0 

t  t  t  t 


subject  to 


Vt 


-  bt  +  Bt-Ixt-1 


Ftxt 


>  f 
—  t 


:  o 


Ltxt  +  eet  *‘t 


=  ft. 


In  these  problems,  b^  represents  the  original  term  b^  from  (P)  as 
modified  by  the  inventory  Bt_^x^  ^  supplied  by  the  current  period  t  -  1 
solution  x<^_j  •  The  dependence  of  b  on  x^  ^  will  usually  be  sup¬ 
pressed  for  clarity.  The  constraints  F  x  >  f  and  Lx  +  e6  >  C 

t  t  —  t  ttt-t 

are  necessary  conditions  for  a  solution  x  to  not  lead  to  future  in¬ 
feasibilities  or  non-optimal  solutions.  Accordingly,  they  are  called 
feasibility  and  look-ahead  (optimality)  cuts,  respectively.  The  constraints 
A  x  =  b  -*■  arc  called  the  body  of  the  constraints. 

The  full  master  problem  (D^)  ,  which  includes  all  possible  proposals 
from  the  future,  is  equivalent  to  (P)  .  Computationally,  a  restricted 
master  is  maintained  at  each  period  by  including  proposals  as  they  are 
generated  and  occasionally  purging  old  proposals  which  are  no  longer  basic  in 
(D  )  .  It  will  be  clear  from  context  which  sets  of  proposals  are  indi¬ 
cated  by  (F^,  f{)  and  (Lt,  f  )  . 
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it  is 


While  nested  decomposition  of  (D)  produces  the  problems  (f)  )  , 
usually  more  instructive  and  more  convenient  computationally  to  work  with 
their  respective  duals  (P  )  .  In  general,  the  problem  at  period  t  acts 

as  a  subproblem  with  respect  to  periods  1 . t  -  1  ,  and  as  a  restricted 

master  problem  with  respect  to  periods  t  +  1,  T  .  The  communication 

between  the  problems  (P  )  can  be  represented  schematically: 


naster  I 


1  \ 

I  ' 

(Vr 

su  h 


\ 


•tits 


I 

modified  j 


\ 


cuts 

modified  \ 

rhs  n~ . . 1 

"xj°w 


Piaster 


Xote  that  the  form  of  the  look-ahead  cuts  is  such  that  they  essentially 
modify  the  objective  function  in  each  (P  )  . 

The  act  of  solving  one  of  the  restricted  master  problems  (PI  is 
called  a  step.  if,  at  some  step,  (P  )  is  infeasible,  a  feasibility 
cut  is  generated  and  imposed  on  (P  )  .  If  a  step  terminates  with  a 
finite  optimal  solution  to  (P  )  ,  a  modified  right-hand  side  for  (Pi 
is  generated.  If  in  the  latter  case  it  is  also  found  that 


z 

t 


>  0 


0 

t-1  ’ 


(l) 


*  0 
where  z ^  is  the  new  optimal  objective  for  (Pt)  and  0^  corresponds 

to  the  most  recent  solution  to  (P^  ,  a  look-ahead  cut  is  generated 

which  may  be  imposed  on  (Pt  .  If  (P  )  is  found  to  have  a  class 
of  solutions  with  objective  unbounded  below,  a  solution  and  a  homo¬ 

geneous  solution  h^  are  generated  in  the  usual  fashion.  The  desired 
effect  of  providing  (P  )  with  a  right-hand  side  of  the  form 


bt+l  +  Vc  +  aBtht 


a  >  0 


(3) 


is  achieved  by  introducing  into  (P  )  a  surrogate  activity,  with  level 
a  _>  0  ,  represented  by  the  column 


"Btht 


(4) 


with  cost  coefficient  c  ht  •  If  (P)  has  a  class  of  solutions  with  ob¬ 
jective  unbounded  below,  eventually  a  ray  indicating  this  will  be  gener¬ 


ated  in  (PT)  .  Otherwise,  a  look-ahead  cut  will  be  generated  in  (pti)  . 
for  some  t'  >  t  ,  which  "cuts  off"  the  ray  successively  in  (PtI  ), 

....  (Pt)  . 


A  wide  variety  of  computational  strategies  may  be  employed  within  the 
framework  described  above.  There  is  freedom  both  in  the  order  in  which 
the  problems  (?t)  are  solved  and  in  when  to  pass  information  between 
problems  in  the  form  of  cuts  and  modified  right-hand  sides.  Computational 
experience  has  indicated  that  the  rate  of  convergence  of  the  algorithm 
can  vary  significantly  when  different  strategies  are  employed.  This  ex¬ 
perience  has  also  indicated  that  in  order  to  attain  the  ease  of  solution 
initially  envisioned  for  this  approach,  signficant  modifications  to  the 
algorithm  described  are  necessary. 
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3 .  A  Modified  Nested  Decomposition  Approach 


3.1  Passing  Surrogate  Columns  Forward 

Suppose,  at  some  step  In  the  course  of  executing  the  algorithm  des¬ 
cribed  in  Section  2.3,  a  finite  optimal  solution  to  (P  )  is  obtained. 

The  optimal  basis  must  include  8^  and,  assuming  A^_  is  of  full  rank, 
at  least  of  the  variables  x  .  Some  slack  variables  corresponding 

to  cuts  which  have  been  imposed  may  also  be  basic.  Let  k^  _>  0  and 
k^  _>  1  be  the  number  of  feasibility  and  look-ahead  cuts,  respectively, 
whose  slack  variables  are  not  basic.  Then  +  k^  +  k^  -  1  of  the 
variables  x  are  in  the  optimal  basis.  Ordinarily,  the  optimal  solution 
x^  is  used  to  form  a  right-hand  side 


bt+l  + 


(5) 


for  the  bodv  of  the  constraints  in  (P  ,)  . 

t+1 

A  modification  to  this  technique  is  outlined  as  follows.  Let 


k  «  +  k2  -  1  (6) 

be  the  number  of  "surplus"  variables,  and  let  0^  be  the  optimal  basis, 
excluding  slacks  from  cuts.  Partition 


6t  =  LBt  U  NBt  U  *  0  } 


(7) 


where  LB 


“t  >  |  NB 1 1  =  k  ,  and  A^g  exists,  where  A^  is  a  con¬ 
venient  abuse  of  the  proper  notation  (Af)  )  •  The  variables  R 

are  called  the  local  basis.  Solve  the  body  of  the  constraints  in  (P  ) 
versus  the  right-hand  side  b  to  yield  the  locally  basic  solution 


-1 
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(8) 


<XlBt)0  *  ALBtEt 


Represent  the  remaining  basic  variables  in 


terms  of  the  local  basis. 


be  the  set  of  homogeneous  solutions  to  the 
are  generated  by  the  representations  Y 

o  c 

locally  basic  solution  (x  )  to  form  a 

LBt 

Vl  +  Bt(xLBt)0 


body  of  the  constraints  which 
.  In  contrast  to  (5),  use  the 
right-hand  side 


(12) 


for  <Pt+1>  • 

gate  activities 


In  addition, 

,  with  levels 


introduce  into  (P  a  set  of  k  surro- 

at+^  0  ,  represented  by  columns  with 


St+1  -  -Vt 


(13) 


in  the  body  of  the  constraints,  zero  coefficients  in  the  cuts  on  (P  ^) , 
and  cost  coefficients  st+^  *  ctHt  '  Pr°blem  is 
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<pt+l> 

minimize  +  c^x^  +  0£+1 


subject  to 


st+lat+l  +  At+lxt+l 

5t+l 

:  71 1+1 

Ft+lxt+l  > 

ft+l 

:0t+l 

l  x  ,  +  e  > 

l 

:  0 

t+i  t+i  t+i  - 

t+1 

pt+l 

“t+r  xt+i  ^  0  • 

essence,  endows  (P  with 

a  right- 

•hand  side 

in  terms  of  the  k  variables 

’ 

and  allows 

nBt 

(Pt+j)  to  select  the  coefficients  at+^  °f  the  parametrization .  Set¬ 
ting  at+^  =  corresponds  to  the  right-hand  side  (5)  obtained  in 

the  pure  nested  decomposition  framework.  Passing  activities  down  the 
staircase  and  allowing  a  ^  to  vary  from  these  values  provides  an 
avenue  of  implicit  communication  between  (P  )  and  (Pj.)  ,  •••»  (P^) 
which  decreases  the  number  of  steps  needed  to  obtain  an  optimal  coordi¬ 
nation  of  the  single-period  problems. 

There  are  two  potential  disadvantages  inherent  in  this  approach. 

First,  it  seems  that  forcing  (P_.,)  to  have  some  or  all  of  the  variables 

t+1 

a  in  its  basis  would  unacceptably  limit  the  number  of  basic  variables 
chosen  from  xt+-^  •  Second,  the  framework  sketched  above  allows  (P 
to  use  the  surrogate  activities  at  any  nonnegative  levels.  In  order  to 
satisfy  the  body  of  the  constraints  in  (P  )  ,  the  relationship 


XI.Bt(cW  =  (xLBt)0-  (YNBt)at+l 


(14) 


must  be  maintained.  Since  (P  )  may  include  surrogate  activities  inherited 
from  (P  ,  in  general  a  parametrization  of  all  variables  which  are  locally 
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basic  in  (P^.) . (P^)  obtained  in  terms  of  the  surrogate  activities 

received  by  (Pt+^)  .  Clearly,  nonnegative  levels  of  at+^  could  cause 

the  parametrized  values  of  variables  in  earlier  periods  to  have  negative 

components.  These  two  points  are  addressed  in  the  next  sections. 

3.2  Degeneracy  in  Optimal  Solutions  Generated 
by  Nested  Decomposition  of  (D) 

It  is  well-known  that  bases  for  (P)  inherit  the  staircase  structure 

of  the  general  problem  and  exhibit  a  ''surplus-shortage"  property  which 

generalizes  the  fact  that  a  basis  which  contains  +  k  of  the  variables 

x^  ,  with  k  >  0  ,  must  exhibit  a  "shortage"  of  k  ,  relative  to  the  num- 

T 

ber  of  remaining  constraints,  i  m  ,  in  the  number  of  basic  variables 

t=2  t 

selected  from  x^,  ...,  x^  .  Fourer  [4]  has  developed  a  set  of  bounds 
on  the  magnitudes  of  the  surpluses  and  shortages  which  each  period  of  a 
basis  may  possess. 

The  following  result  describes  a  manifestation  of  this  property  in  the 
setting  of  the  algorithm  developed  in  Section  2.3. 

Theorem.  Suppose  an  optimal  coordination  of  the  single-period  problems 
(P^)  has  been  obtained.  If  t  <  T  and  (Pt)  has  k  +  1  cuts  whose 
slack  variables  are  nonbasic,  then  the  basic  solution  to  (P^+^)  has  at 
least  k  degenerate  variables. 

Under  the  modification  outlined  in  Section  3.1,  surrogate  columns 
passed  from  period  t  replace  degenerate  variables  in  the  optimal  basis 
in  period  t  +  1  .  These  degenerate  variables  may  include  basic  slack 
variables  for  the  cuts  in  period  t  +  1  .  The  fact  that  the  surrogate 
columns  replace  degenerate  variables  guarantees  that  no  solutions  of 
interest  are  excluded  in  (Pt+^)  . 
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3 . 3  Maintaining  Feasibility  of  Locally 
Basic  Solutions 

As  indicated  at  the  end  of  Section  3.1,  nonnegativity  of  does 

not  imply  nonnegativity  of  the  parametrized  values  of  the  variables  which 
are  locally  basic  in  periods  1,  t  .  Two  general  approaches  to  sur¬ 

mounting  this  difficulty  are  possible.  First,  the  feasibility  of  earlier 
periods  can  be  ignored  during  the  optimization  of  the  modified  problem 
(P  and  restored  in  a  subsequent  procedure  which  would  minimize  an 

appropriate  infeasibility  form.  The  alternative  is  to  employ,  while  op¬ 
timizing  (P  j)  ,  an  extended  minimum-ratio  test  which  follows  the  para¬ 
metric  variation  of  locally  basic  variables  in  earlier  periods  and  indi¬ 
cates  when  such  a  variable  blocks  the  increase  of  an  incoming  column  in 

The  latter  approach  is  adopted  here.  The  representation  of  an  incoming 
column  in  (P  includes  weights  on  any  surrogate  columns  which  are  in 

the  local  basis.  Using  these  weights  and  repeatedly  applying  (14)  and  (9) 
yields  a  representation  of  the  incoming  column  in  terms  of  the  variables 
which  are  locally  basic  in  periods  1,  ...,  t  +  1  .  This  representation 
is  used  to  implement  the  extended  minimum-ratio  test.  This  scheme  of 
local  inverses  linked  by  representations  of  surrogate  columns  can  be 
viewed  as  a  compact  basis-inverse  technique  which  maintains  a  nearly 
block-angular  inverse  of  the  columns  which  are  locally  basic  in  periods 
1,  ....  t  +  1  . 

When  the  extended  minimum-ratio  test  reveals  that  a  variable  in  a 
period  prior  to  t  +  1  blocks  the  increase  of  an  incoming  column  in 
(P  several  strategies  may  be  employed.  The  choice  in  this  work  is 
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to  "shuffle"  the  structure  of  the  surrogate  columns  by  exchanging  the 
roles  of  some  of  the  activities  which  are  and  are  not  locally  basic. 

This  process  is  designed  to  take  the  blocking  variable  and  pass  it  down 
to  (Pi)  as  a  surrogate  column.  Primal  and  dual  solutions  in  each 
period  are  unchanged,  and  when  optimization  of  (P  is  resumed,  the 

same  incoming  variable  is  blocked  by  this  surrogate  column.  The  indi¬ 
cated  pivot  in  (P  is  made  and  maintains  the  nonnegativity  of  locally 

basic  solutions  in  periods  prior  to  t  +  1  . 

3 . 4  Computational  Strategies  and  Further  Work 

Again,  a  wide  variety  of  computational  strategies  may  be  employed 
within  the  framework  described  above.  There  is  freedom  in  the  order 
in  which  the  problem  (P^)  are  solved  and  in  when  to  pass  information 
between  problems  in  the  form  of  cuts,  modified  right-hand  sides,  and 
surrogate  columns.  The  presence  of  the  surrogate  columns  opens  addi¬ 
tional  options,  including  variations  of  the  "shuffle"  described  in  Sec¬ 
tion  3.3. 

Currently,  computational  experience  is  being  obtained  with  a  code 
written  by  Wittrock  in  Mathematical  Programming  Language  at  Stanford 
University.  This  language  facilitates  experimentation  of  the  type  needed 
at  this  stage  of  the  work.  The  thrust  of  this  experimentation  is  to  de¬ 
vise  computational  strategies  which  tend  to  minimize  the  computational 
effort  needed  to  obtain  an  optimal  solution  to  (P)  .  Since  the  manipu¬ 
lations  of  data  structures,  and  the  form  and  frequency  of  updates  to  the 
local  inverses,  depend  heavily  upon  the  computational  strategies  which 
are  employed,  decisions  about  these  factors  have  not  yet  been  made. 
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